Supermagnetosonic Jets behind a Collisionless Quasi-parallel Shock 
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The downstream region of a collisionless quasi-parallel shock is structured containing bulk flows 
with high kinetic energy density from a previously unidentified source. We present Cluster multi- 
spacecraft measurements of this type of supermagnetosonic jet as well as of a weak secondary shock 
front within the sheath, that allow us to propose the following generation mechanism for the jets: The 
local curvature variations inherent to quasi-parallel shocks can create fast, deflected jets accompanied 
by density variations in the downstream region. If the speed of the jet is super (magneto) sonic in 
the reference frame of the obstacle, a second shock front forms in the sheath closer to the obstacle. 
Our results can be applied to collisionless quasi-parallel shocks in many plasma environments. 

PACS numbers: 52.35.Tc, 52.40.Kh, 96.50.Ek 



Introduction. — When the angle between the nominal 
shock normal and the upstream magnetic field is small, 
the shock transition in a collisionless plasma is much 
more complex than in the quasi-perpendicular case []]]. 
The nonthermal nature of the upstream side of a quasi- 
parallel shock has been recognized for decades @, [3, El • 
The downstream region, however, has only recently come 
under active research, both in astrophysical (supernovae 
[H) and Solar System (termination shock 0, Earth's bow 
shock 0,0]) contexts. 

The most detailed and extensive data of collisionless 
shock waves are from the Earth's bow shock. In con- 
trast to remote observations and laboratory measure- 
ments, the near-Earth space can be used to study in 
situ supersonic plasma flow past a magnetic obstacle — 
the flow of the solar wind around the magnetosphere of 
the Earth. The magnetospheric boundary (the magne- 
topause) is usually located at a distance of 10 Earth radii 
(1 i?E = 6371 km) in the solar direction. The bow shock 
is curved at magnetospheric scales while the structures in 
the solar wind and interplanetary magnetic field are large 
compared to the size of the magnetosphere. Hence the 
locations of parallel and perpendicular regions of the bow 
shock vary depending on the direction of the interplane- 
tary magnetic field. Consequently, we can access a wide 
range of plasma conditions via spacecraft observations. 

Recent measurements have shown that the flow in the 
downstream region of a quasi-parallel shock is structured: 
Nemecek et at [9J] have reported observations of transient 
ion flux enhancements in the Earth's magnetosheath dur- 
ing radial interplanetary magnetic field. In subsequent 
studies, Savin et al. [Io| have found more than 140 events 
of anomalously high energy density. However, the source 
of these jets of high kinetic energy and ion flux has re- 



mained unclear. In this Letter, we present a set of multi- 
spacecraft measurements from Cluster [ll| that allows us 
to suggest a formation mechanism for such jets. 

Data. — We have analyzed Cluster measurements from 
the evening of March 17, 2007, when the four spacecraft 
(C1-C4) were close to the nose of the magnetosphere. 
The spacecraft constellation was quite flat in the nom- 
inal plane of the magnetopause, since C3 and C4 were 
close to each other (950 km, 0.15 Re apart), while the 
others were slightly more than 7000 km (>1 Re) away. 
We have used data from the magnetic field experiment 
FGM from all four spacecraft, and from the ion experi- 
ment CIS-HIA from CI and C3 [11| . Information about 
the upstream conditions was provided by ACE and Wind 
satellites situated near the Lagrangian point LI, as well 
as the Geotail spacecraft, which at the time was in the 
foreshock region near the subsolar point. 

The free upstream solar wind flow was quite fast 
(V ~ 530 km/s) and steady (see the upper panels 
of Figure [1]) . The particle number density was around 
2 cm~ 3 and hence the dynamic pressure (pV 2 , where 
p is the mass density) was low, close to 1 nPa. The 
interplanetary magnetic field was approximately radial, 
i.e., the sunward magnetic field component Bx [22| was 
dominant. Moreover, the angle between the flow direc- 
tion and the magnetic field was less than 20°. Conse- 
quently, the bow shock was quasi-parallel at the subso- 
lar point. The upstream Mach numbers [23| were all 
larger than 10: M A ~ 12, M s ~ 16, and M M s ~ 10. 
The location of the bow shock, as observed by Geotail at 
(X,Y 1 Z) GSE = -(14, -7, -3) R E [22] when the shock 
moved over the spacecraft several times between 17:30 
and 24:00 UT, matches well to the empirical model [12| 
for the measured upstream parameters. 
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FIG. 1: Upper panels: upstream solar wind data from the 
ACE satellite (time-shifted by 44 minutes to account for the 
solar wind propagation to the magnetopause) . First panel: 
magnitude of the interplanetary magnetic field and angle 9 
between the field direction vector and the XGSB-axis [22]. 
Second panel: plasma number density and dynamic pressure. 
ACE was located at (X,Y,Z) GS e = (237, 36.4, -18.6) R E . 
The gray shading marks the period of interest. Lower panels: 
magnetic field from all four Cluster spacecraft in GSE coor- 
dinates. The quartet was situated around (10.7, 1.5, 3) -Re- 
The color panels mark different plasma regions. White back- 
ground between color panels represents transition between 
two regions. 



The Cluster quartet, moving on an outbound or- 
bit near the subsolar point, encountered the magne- 
topause the first time shortly after 17:00 UT and passed 
into the solar wind at 22:30 UT. Between 17:00 and 
20:00 UT Cluster observed multiple magnetopause cross- 
ings. Moreover, during this 3-hour period Cluster ob- 
served several high speed jets (V ~ 500 km/s) in the mag- 
netosheath behind the quasi-parallel bow shock. Here we 
concentrate on the jet between 18:13 and 18:17 UT. 

As displayed in the lower panels of Figure [TJ all four 
spacecraft were inside the magnetosphere at the begin- 
ning of the interval. First the magnetopause moved in- 
wards passing over the Cluster quartet at 250 km/s (ob- 
tained using 4-spacecraft timing). Then the spacecraft 
observed a weak shock within the magnetosheath mov- 
ing in the same direction at 140 km/s. In the first panel 
of Figure [21 the CI measurements show that at this mo- 
ment the component of the plasma velocity parallel to 



the shock normal in the reference frame moving with the 
shock V n exceeds the magnetosonic speed Ums (as well 
as the other characteristic speeds [23]). Hence the mag- 
netosonic Mach number Mms > 1- (This is also the case 
with respect to the magnetopause.) The same was ob- 
served by C3 located about 8000 km away (not shown). 

After the shock Cluster entered a cold, supermagne- 
tosonic jet with a plasma speed close to 500 km/s (see 
Figure^ first panel). At the location of C2, the shock 
and the magnetopause moved back across the spacecraft 
and it re-entered the magnetosphere for several seconds 
at 18:15:20 UT. The other spacecraft stayed in the su- 
personic jet for over a minute moving gradually back into 
normal sheath-type plasma. This transition can be seen 
in the ion velocity distributions (not shown) : the narrow 
(~1 MK) distribution of the jet was slowly replaced by 
a warmer, symmetric quasi-Maxwellian after 18:16 UT. 

While in the jet, Cluster observed a gradual increase in 
both plasma density and magnetic field: from low values 
of 7 cm~ 3 and 8 nT at the beginning of the jet to very 
high values of 22 cm -3 and 30 nT at the end (see Fig- 
ure [2 second panel for CI ion density). Consequently, 
the dynamic pressure in the jet increased to over 6 nPa, 
as compared to the nominal pressure of 1 nPa. These en- 
hancements were accompanied by a substantial deflection 
of the bulk flow from its nominal direction, as illustrated 
by the third panel of Figure El 

Interpretation. — We propose the following mechanism 
to explain the formation of the jet: First, consider an 
oblique shock with radial upstream conditions (Vi || Bi) 
as illustrated in the inset of Figure [3] The Rankinc- 
Hugoniot jump conditions for high Ma give V\ n = r~V<i n 
and Vu ~ Vzt, where r is the shock compression ratio. 

We then consider the streamlines of plasma flow across 
a curved high Ma shock as illustrated in Figure [3] We 
infer, based on the analysis of the observations as will be 
discussed below, that the scale of the shock ripple under 
consideration is of the order of the spacecraft separation: 
50 - 100 ion inertial lengths, 7000 - 15000 km, 1 - 3 R E . 
As the shock primarily decelerates the component of the 
upstream velocity Vi parallel to the shock normal n, the 
shock crossing leads to efficient compression and deceler- 
ation in regions where the angle a between Vi and n is 
small. Wherever a is large, however, the shock mainly 
deflects the flow while the plasma speed stays close to the 
upstream value. The plasma is still compressed so that 
the higher density together with the high speed leads to 
a jet of very high dynamic pressure. Furthermore, if the 
speed V% of this jet on the downstream side is still su- 
per(magneto)sonic in the reference frame of the obstacle, 
a second shock front forms closer to the obstacle. In addi- 
tion, depending on the ripple geometry, the flow behind 
the shock can converge causing local density enhance- 
ments, or diverge causing density depletions. 

Let us compare Figure [3] with the CI measurements 
presented in Figure [2] where, in the third panel, the bulk 
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FIG. 2: First panel: Component of plasma velocity paral- 
lel to the local shock normal (—0.59, 0.52, — 0.61)gse of the 
secondary shock (calculated with minimum variance analy- 
sis), in the reference frame moving with the secondary shock 
V n (black solid curve) and total plasma speed V (light blue 
solid curve). Dashed curves show the characteristic speeds: 
Alfven speed Va (blue), sound speed Vs (green), and magne- 
tosonic speed Vms (red) [23]. Second panel: Plasma number 
density. Third panel: Bulk velocity projection to ( — Zgse, 
-^GSe) plane. Fourth panel: The angle a calculated from 
the observed velocity deflection using the jump conditions for 
high A/a and r = 4. The calculation is not expected to be 
valid at the edges of the jet where the shock is weak, and 
hence a is shown for the center only. All data are from CI. 
The color coding for different plasma regions is the same as 
in Figure [TJ 
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FIG. 3: Illustration of the effect of local shock curvature. The 
variation of the plasma number density in the downstream re- 
gion is illustrated by the shading: dark blue indicates density 
enhancement, light blue indicates density depletion. The tra- 
jectory of CI in the reference frame moving with the ripple 
is sketched with the dashed line. The inset details the flow 
deflection when Vi is not parallel to n. 



flow direction is displayed in the (— Zqse, -^gse) pl ane - 
The observed pattern of the supermagnetosonic flow af- 
ter the secondary shock suggests that there is a ripple 
in the bow shock similar to the one of the illustration 
moving in the ~ Zqse direction. This interpretation is 
supported by the observed density and flow speed pro- 
files. The fourth panel of Figure [2] shows the upstream 
angle a for the supermagnetosonic jet calculated from the 
observations (considering both the downstream and up- 
stream data and taking r = 4). During the main velocity 
deflection, a ~ —65°. The flow pattern in the Ygse (not 
shown) reveals more of the three-dimensional structure 
of the ripple and will be considered elsewhere. The ob- 
servations of C3 are similar, though not identical to CI. 
Given this and the fact that C2 was outside of the jet, we 
infer that the lower limit for the scale of the bow shock 
perturbation is of the order of the spacecraft separation 
(~ 8000 km, 1.2 Re). 

The ripples we propose to be the source of the high 
speed jets stem from the unstable nature of collisionless 
quasi-parallel shocks: reflected ions can stream against 
the upstream flow and interact with the incident plasma 
over long distances before returning (if at all) to the 
shock. This interaction triggers instabilities and creates 
waves that steepen into large structures convecting back 
to the shock front (see [l[, and the references therein). 
The effect is most pronounced when Bi and Vi are 
aligned in the coordinate system of the obstacle. 

Both observations and simulations have shown that 
ripples are inherent to quasi-parallel shocks: Observa- 
tions of the ion reflection on the upstream side of the 
Earth's bow shock Q indicated that the direction of n 
varies when the shock is quasi-parallel at the subsolar 
point. Such studies on the ion distributions have also 
shown that, at times, the solar wind does indeed pass 
through the shock layer without significant heating [13j. 
However, no connection between these two findings was 
made. Furthermore, recent multi-spacecraft observations 
have characterised in detail the short, large amplitude 
magnetic structures (SLAMS) [3, EH convecting in the 
upstream of Earth's quasi-parallel bow shock towards the 
shock front. SLAMS have a scale size up to 1 Re compa- 
rable to the ripples discussed here. In addition, measure- 
ments showed signatures that the shock transition itself 
is narrow consisting of only one to a few SLAMS. The 
roughness of the parallel part of the shock front due to 
SLAMS is clearly seen in the bow shock simulations of, 
e.g., Blanco-Cano et al. [161 ]. 

Discussion and Conclusions. — In previous studies 
Savin et al. [l(| have found several jets with high ki- 
netic energy density (^pV 2 ), of which 33 jets had an 
energy density larger than 10 keV/cm 3 (compared with 
19 keV/cm 3 in this event). Nemecek et al. [9] have 
also reported what they call transient ion flux enhance- 
ments, with fluxes of 6 x 10 8 /cm 2 s (~ 8 x 10 8 /cm 2 s in 
this event), during intervals of radial interplanetary mag- 
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netic field. Both transient flux enhancements and high 
kinetic energy jets have properties similar to the jets re- 
ported here. Neither Nemecek et al. nor Savin et al. 
could identify a clear source for the jets, but they could 
rule out, e.g., reconnection. Here, we propose a gener- 
ation mechanism for high speed jets in the sheath that 
is in agreement with the measurements presented in this 
Letter and those reported in previous studies. Naturally, 
we cannot ascertain that all of the previously reported 
jets stem from the same shock geometry related origin. 

It has become evident that shocks are more structured 
than was previously recognized, so that a conventional 
plane wave description is not sufficient. In fact, the mech- 
anism proposed in this Letter for spatial structuring of 
the downstream is valid for all rippled shocks regard- 
less of magnetic field obliquity, provided that the Mach 
number is high. As Voyager 1 and 2 crossed the helio- 
spheric termination shock [6[ , their observations revealed 
a rippled, supercritical (Mms ~ 10) quasi-perpendicular 
shock [Tfl] . Likewise, interplanetary shocks seem to be 
nonplanar [T3| and also oblique ones may be rippling [h| . 
Therefore we expect that the effects of the ripples, includ- 
ing supersonic jets, can be observed behind collisionless 
shocks in many plasma environments, and especially be- 
hind extended, varying shock fronts having quasi-parallel 
regions. 

In astrophysical context, the high speed jets and 
nonthermal structure can act as seeds for magnetic 
field amplification and particle acceleration [Hj], even 
for smooth upstream plasma. In magnetospheric con- 
text, the jets with their high dynamic pressure pro- 
vide a previously unidentified source for magnetopause 
waves during steady solar wind conditions. A locally 
perturbed magnetopause is consistent with the Cluster 
measurements of C2 being within the magnetosphere 
while the other spacecraft were in the jet. In turn, the 
large magnetopause perturbation can affect the coupled 
magnetosphere- ionosphere dynamics [2(| • Note also that 
this Letter presents observations of a weak shock within 
the magnetosheath during steady upstream conditions. 
Previous studies of discontinuities within the sheath have 
been related to bow shock interaction with interplanetary 
shocks (see 21 1, and the references therein). 

In summary, we propose a generation mechanism for 
high speed jets in the downstream side of a quasi- 
parallel shock based on a set of multi-point measure- 
ments. Quasi-parallel shocks are known to be rippled 
even during steady upstream conditions. The local cur- 
vature changes of the quasi-parallel shock can create fast 
bulk flows: in the regions where the upstream velocity is 
quasi-perpendicular to the local shock normal, the shock 
mainly deflects plasma flow while the speed stays close 
to the upstream value. Together with the compression of 
the plasma, these localized streams can lead to jets with 
a kinetic energy density that is several times higher than 



the kinetic energy density in the upstream. 
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